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Abstract. Quantum effects in statistical mechanics are important when the thermal 
wavelength is of the order of, or greater than, the mean interatomic spacing. This is 
examined at depth taking the example of a hard-sphere Boltzmann gas. Using the 
virial expansion for the equation of state, it is shown that the interaction energy of a 
classical hard-sphere gas is exactly zero. When the (second) virial coefficient of such 
a gas is obtained quantum mechanically, however, the quantum contribution to the 
interaction energy is shown to be substantial. The importance of the semiclassical 
corrections to the interaction energy shows up dramatically in such a system. 
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1. Introduction 

The virial expansion of a dilute gas is standard fare in the curriculum of the statistical 
mechanics course for beginning graduate students. The hard-sphere classical gas 
provides a nice example on this topic, since the first few classical virial coefficients are 
analytically known. The quantum corrections in the dilute hard-sphere gas, however, is 
a topic not covered in text books. 

The objective of this paper is to study the quantum corrections in a hard-sphere 
one-component Boltzmann gas. The usefulness of semiclassical approximations in 
understanding the difference between the classical and quantum results is demonstrated 
in this paper. As usual, we consider the gas in the thermodynamic limit, where the 
number of atoms (of one kind only), N — > oo, as also the volume of the enclosure, 
V — > oo, but the number density n = N/V remains finite. Each atom is taken to 
be a rigid (impenetrable) sphere of diameter a, so that the minimum relative distance 
between the spheres is a. It must be mentioned at the outset that this system has been 
studied extensively in the last fifty years (see [T] for a nice summary). Semiclassical 
quantum corrections in the hard-sphere Boltzmann gas were examined (amongst others) 
by Hill [2j and Gibson 3J. Our focus in this paper will be mainly on the interaction 
energy of the hard-sphere gas. We first show, using the virial expansion for the equation 
of state, that this is exactly zero when the virial coefficients for the hard-sphere gas 
are calculated classically. Since we consider a Boltzmann gas, there are no exchange 
contributions due to statistics. For a dilute gas at high temperature, it is the second 
virial coefficient that gives the dominant contribution. The second virial coefficient 
may also be calculated quantum mechanically from the Beth-Uhlenbeck formula using 
the hard-sphere scattering phase shifts in the various partial waves [3]. Contrary to 
naive expectations, the quantum result for the interaction energy of the Boltzmann 
gas at high temperatures is substantially different from the classical one, even when 
the thermal wavelength is only a fraction of the hard-sphere diameter. Obtaining the 
quantum result requires careful numerical work, adding contributions from more and 
more partial waves with increasing temperature. The difference between the quantum 
and the classical results may be accounted for by the semiclassical corrections, given in 
the form of an asymptotic series. This example shows the usefulness of the semiclassical 
method in statistical mechanics even for a singular interaction. 

2. The interaction energy 

2.1. The virial expansion 

To appreciate this in more detail, note that the virial expansion of a (one-component) 
gas (classical as well as quantal) for the pressure P at a temperature T is given by 




(1) 
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where 02,03, . . . , are the dimensionless second, third, etc., virial coefficients, r = /c#T, 
n = N/V is the number density of particles, and A = ^2irh 2 /mT is the thermal 
wavelength. Even for an ideal quantum gas of bosons or fermions, the virial coefficients 
are nonzero. These contributions due to statistics will not be included in our Boltzmann 
gas; only the part due to the interaction between the atoms will be considered. The free 
energy F = E — TS may be easily obtained by integrating the above P with respect 
to the volume V (since P = —(dF/dV) T ), the constant of integration being chosen to 
match the perfect Boltzmann gas value: 
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The energy may next be obtained by using the relation E = —T 2 [d(F/T)/dr}. After 
subtracting out the energy of the perfect gas part, one obtains the virial series for the 
interaction energy per particle: 
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The above equation for the interaction energy, modified for a two-component Fermi gas, 
has been used recently by Ho and Mueller jH]. Whereas for a classical gas the virial 
coefficients may be expressed as integrals involving the interaction potential, for the 
quantum problem a solution of the j-body problem is needed to obtain aj. Henceforth 
we shall denote the classical virial coefficients by Aj to differentiate from their quantum 
counter parts, to be still denoted as aj. 

It is easy to see now why the interaction energy for a classical hard-sphere gas is zero 
order by order, in every power of (nX 3 ). For a hard-sphere diameter a, A 2 = ^(a/A) 3 
(see (jnj) below), and the higher order Aj's may be expressed as Aj = K(A 2 y~ 1 , where 
n is a constant 1J. Substituting this in (jSJ), we see that the interaction energy of 
the classical hard-sphere gas vanishes identically in each power of (nX 3 ). In fact, the 
classical hard-sphere gas is rather special, since the right hand side of (0) is independent 
of temperature. 

We note that in the regime (nX 3 ) « 1, the series (j3J) may be terminated after the 
first order term 
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The classical second virial coefficient A 2 is given by [3] 
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where (3 — 1/r, and V(r) is the interatomic two-body potential. The quantum second 
virial coefficient requires a knowledge of the bound-state and the continuum spectra. 
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For a one-component quantum gas, 6 2 = — a 2 is given by jSj 
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In the above, E n /s are the binding energies of the two-body bound states (absent in 
our problem), and Si is the scattering phase shift in the Zth partial wave. It is important 
to note that our (@) differs from the usual definition given for quantum gases in text 
books 0[7|: the overall multiplicative factor in our case being a/2 rather than 2y/2. 
The second virial coefficient, after all, is just the two-body partition function for relative 
motion, and for the Boltzmann case it has a factor 1/N\, which explains why our factor 
is less by 2 for the two-body case. Our definition here is in agreement with the earlier 
authors El E] who considered the Boltzmann case. 

2.2. Numerical comparison 

We compare the classical second virial coefficient A2 = ^-(a/A) 3 with its quantum 
counter part — 6 2 as defined by © for the hard-sphere gas. The phase shift 5i for 
hard-sphere scattering is given by [Oj tan 5/ = ji(ka) / ni(ka) , and a little algebra shows 
that 

d5i _ __a 1 

dk (ka) 2 jf(ka) +nf(ka) 

For a fixed a, as one goes to higher temperatures (i.e., larger values of (a/A)), more and 
more / values are needed to get convergence in the computed b$ using (jBJ). For example, 
for (a/A) = 2, convergence is obtained only when the sum extends to Z max = 33 ^f. This 
may be understood by noting that for small (3 (high temperature), the damping of 
the integrand from the gaussian is less effective, and higher values of the variable k 
are needed for the convergence. From impact parameter argument ^U], the cut off 
Z max ~ (ka) in the sum over I therefore also increases. In figure ^ we display the 
classical and quantum second virial coefficients (A2 and —62) as a function of (a/ A). 
One expects the classical regime for (a/ A) » 1. As is clear from figure ^ however, 
the two remain very different even in the so called classical regime. This difference may 
be understood from semiclassical corrections, as we shall see in the next section. For a 
perfect Boltzmann gas, the energy per particle, E /N = |r. In figure El we plot E int /E 
in units of (nA 3 ). Recall that this this quantity is exactly zero classically for any value 
of (nA 3 ). In figureEl the interaction energy is calculated from (J3J) to first order in (nA 3 ), 
using the virial coefficient (—62) given by ©• The disagreement between the classical 
and the quantum values of the second virial coefficient, as also the interaction energy, 
keeps increasing for larger values of (a/ A). This, however, is misleading, as we see below. 

% The quantity hi was calculated by first summing over I and then integrating. We specified a minimum 
relative error in the integrand (or the sum) of 1CP 3 . 
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Figure 1. The quantum and classical second virial coefficient with semiclassical 
corrections (to be discussed later.) The lowest curve gives the classical second virial 
coefficient and semiclassical approximations approach the quantum mechanical curve 
for sufficiently large a/ A. 




Figure 2. The quantum and semi-classical approximations to the interaction energy 
in units of n\ 3 . 
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For fixed a and n, larger values of (a/A) imply smaller (nA 3 ), which actually decreases 
the pressure P/(nr) in (JTJ) and the energy ratio E int /E in magnitude. 

To get a better feel for figure note that the maximum density for close packing 
of the hard spheres [T] is n max = -y/2/a 3 . In a dilute gas, let us assume that 
n = n max /10 = V2/(10a 3 ). We then get 

.XJ 10(nA 3 ) ' y ' 

If we assume (nA 3 ) = 0.1, we get (a/A) 3 = a/2, or a/A ~ 1. From figure El we see 
that E int /E ~ 0.12 in this case as opposed to zero if this was calculated classically. 
Semiclassical corrections, to be discussed later, are still important in this regime. On 
the other hand, for a gas at equilibrium at room temperature and pressure, nA 3 is at 
least ten times less, so quantum corrections are negligible. 

In figure El we similarly plot C int /C in units of (nA 3 ) as a function of (aA). 
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Figure 3. The specific heat calculated quantum mechanically and using semi-classical 
approximations. Since the second term in the semi-classical expansion is zero, the one- 
term and two-term approximations are identical. 

Classically, of course, Co is 3/2 per particle in units of the Boltzmann constant, and 



3. The second virial coefficient 

3.1. The classical- quantum connection 

In this section, the connection between the classical .4.2 and its quantum counterpart 
— & 2 (/?) w iH be explored for a given two-body potential V(r). This will help us to 
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understand why the classical and quantum results are so different for the example 
that we are studying. For this purpose, it is convenient to make the partial wave 
decomposition 
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it may be verified easily (by replacing the sum over I by an integral), that 
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One may enquire about the relation between this expression (jlOj) and the quantum 
counterpart ■ This question was actually answered nearly seventy years back by 
Kahn ^1] in his doctoral dissertation. We shall give here a new compact derivation of 
his result, considering only the part of the virial coefficient arising from the scattering 
channel in the continuum. Using (JHJ), and the relation E = h 2 k 2 /m, we have (with no 
bound states) 
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We note that 4=&2^(/3) is j us t the Laplace transform of the derivative of the phase shift 
with respect to the energy. In the lowest order WKB approximation, the phase shift is 
given by [IT 
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In the above, the effective potential VJ(r) is defined as 
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In fll2j) . r\ and r$ are the classical turning points where the respective integrands go to 
zero. The derivative of the phase shift may be written as 
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in terms of the unit step function Q(x). By noting that the Laplace transform of 
— ^1 i s /^"exp^—^^ we immediately obtain the desired result 
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If the lowest order WKB approximation (|14|) for a given potential were exact, the 
classical and quantum (partial wave) virial coefficients would be the same. In general, 
of course, there are higher order semiclassical corrections in powers of % to the WKB 
lowest order term that accounts for the deviation from the classical result. We next 
examine the semiclassical corrections. 



3.2. Semiclassical corrections to A2{0) 



There is the Wigner-Kirkwood ^21 E3 El semiclassical expansion of the quantum 
partition function (or the second virial coefficient) in powers of h 2 and the derivatives 
of the potential. This, however, is not applicable for the hard-sphere potential, since 
its derivatives do not exist. For the hard-sphere problem Hill |2j and Gibson [3 j used 
advanced techniques in scattering theory to obtain the following asymptotic expansion 
for the virial coefficient : 
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where the classical virial coefficient A% = ^(a/A) 3 . Since A is proportional to H, the 
above series is in powers of h, and not h 2 . This is characteristic of potentials with sharp 
surfaces jSj. The derivation of the above series is quite involved, and will not be given 
here. We can present, however, a qualitative argument to indicate why one may expect 
powers of (A/a) in the semiclassical series above. Consider the integral (jHJ for A2{(3) for 
the hard spheres. The potential V(r) is zero whenever the relative distance r between 
the centers of the spheres is greater than a. For r = a the spheres touch each other with 
their sharp surfaces. Semiclassically, though, we may expect a fuzziness (or uncertainty) 
in the contact distance by an amount which is of the order of the de Broglie wave length 
Ad = h/(p), where p is the absolute value of the relative momentum between the two 
spheres. The average (p) with Boltzmann distribution, and the corresponding de Broglie 
wavelength are given by 
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In evaluating the integral (0) semiclassically, we may put V(r) = for r > (a + Xd), 
so that now the second virial coefficient becomes ^(a + A^) 3 /A 3 . This may be written 
in the form of a series in powers of (A/a), multiplied by the classical term. This series, 
unlike the one given by (J16)) . terminates with the cubic term. Although the coefficient 
of the leading (A/a) term is correctly reproduced, the agreement is fortuitous. The 
main merit of this derivation is simply to show that it is not unreasonable to expect 
semiclassical corrections in powers of (A/a). 
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The order by order semiclassical corrections using (J 16)) are computed and shown 
for the second virial coefficient (figure HJ), the ratio of E int /E (figure 0), and the specific 
heat ratio C; n t/Co (figure EJ). We see from these that the semiclassical corrections are 
able to reproduce the quantum behaviour very well for (a/A) > 0.2. Since the series 
fTK|) is an asymptotic one, the higher order terms are seen to diverge for small values 
of (a/A). Nonetheless, the series is very accurate for large values of (a/A). For a given 
(a/A), the asymptotic series under consideration has to be terminated after an optimum 
number of terms are summed. 

We have thus demonstrated the usefulness of the semiclassical series for 
understanding the classical to quantum behaviour. 
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